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We present ab initio calculations of transport properties of atomic-sized aluminum contacts in the 
presence of oxygen. The experimental situation is modeled by considering a single oxygen atom (O) 
or one of the molecules O2 and O3 bridging the gap between electrodes forming ideal, atomically 
sharp pyramids. The transport characteristics are computed for these geometries with increasing 
distances between the leads, simulating the opening of a break junction. To facilitate comparison 
with experiments further, the vibrational modes of the oxygen connected to the electrodes are 
studied. It is found that in the contact regime the change of transport properties due to the presence 
of oxygen is strong and should be detectable in experiments. All three types of oxygen exhibit a 
comparable behavior in their vibrational frequencies and conductances, which are well below the 
conductance of pure aluminum atomic contacts. The conductance decreases for an increasing number 
of oxygen atoms. In the tunneling regime the conductance decays exponentially with distance and 
the decay length depends on whether or not oxygen is present in the junction. This fact may provide 
a way to identify the presence of a gas molecule in metallic atomic contacts. 

PACS numbers: 72.10.-d,71.15.-m,63.22.+m,73.23.-b 



I. INTRODUCTION 

Since the first experiments on atomic contacts, the "at- 
mospheric" surroundings are known to have a strong ef- 
fect on their charge transport properties.^ Thus, studying 
the influence of adsorbed gas molecules is of great impor- 
tance for the design of functional molecular-scale elec- 
tronic components. It also shows some promise for the 
future development of gas sensors for "electronic noses" . 
However, the systematic experimental investigation of 
the influence of single gas molecules on charge trans- 
port through metallic atomic contacts started only a few 
years ago with hydrogen in atomic platinum contacts 
It was found that the zero-bias conductance through 
such platinum-hydrogen contacts tends to be close to the 
quantum of conductance Go = 2e^//i. In addition, this 
conductance value is due to a single, almost fully trans- 
parent channel^^i^ Theoretical work has mostly confirmed 
these findings j^i^i^iii^ 

Even so, one of the main problems in experiments of 
this type is to know if and in which configuration gas 
molecules are present in the contact region. A definite 
identification of the presence and the type of gas can 
be based on statistical techniques, such as the measure- 
ment of conductance histograms For individual con- 
tacts, this becomes much harder. In the experiments of 
Refs. [iljsl the identification of a single hydrogen molecule 
was based on the so-called inelastic point contact spec- 
troscopy, where the local vibrational modes can be seen 
in the current- voltage (/-V) characteristics as abrupt fea- 
tures at bias voltages corresponding to the energies of the 
vibrations.'^ The evolution of these features with stretch- 



ing of the contacts can then be compared to calcula- 
tions of the vibrational modes for given configurations 
of the molecule. In fact, theoretical methods have re- 
cently been developed that enable a direct comparison 
between the full experimental and theoretical I-V char- 
acteristics, including the inelastic signatures of the vibra- 
tional modes^^^^^^ 

For this purpose, namely to identify small molecules in 
atomic metal contacts through transport measurements, 
we study the general properties of oxygen (O) in atomic 
aluminum (Al) contacts. In this paper we concentrate 
on the behavior of the zero-bias transport characteristics 
when the atomic contacts are stretched. In the spirit of 
Refs. [iljslj^, we also study the evolution and character of 
the vibrational modes. Aluminum was chosen, because 
it becomes superconducting below temperatures of 1 K. 
As a result, in addition to measurements of the total 
zero-bias conductance, it is possible to determine exper- 
imentally the transmission eigenvalues of the individual 
channels from the I-V characteristics in the supercon- 
ducting state. As particular examples, we study single 
oxygen atoms (O), oxygen dimers (O2), and ozone (O3). 
These molecules (we also refer to the oxygen atom as a 
molecule) are placed initially in the middle of the con- 
tact between two aluminum electrodes, since their effect 
on the transport properties is then at its largest. We de- 
scribe our contacts at the level of density-functional the- 
ory (DFT). In particular, the DFT electronic structure 
is employed to determine the charge transport properties 
by means of Green's function techniques. 

In the following we present stable structures of the 
three types of contacts, where oxygen molecules bridge 
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the gap between two idealized aluminum electrodes. For 
them, we calculate the vibrational modes and the low- 
temperature transport properties, including the charac- 
terization of the individual transmission channels. We 
find that both the energies of the vibrational modes, as 
well as the evolution of the conductance with increasing 
electrode distance show qualitatively comparable charac- 
teristics for the three investigated molecules. In contrast 
to pure Al, the conductances of the Al-O^^-Al contacts 
are always well below IGq, and for an increasing number 
of O atoms, the conductance decreases. Moreover, the 
exponential decay lengths in the tunneling regime are 
different in the presence and in the absence of oxygen. 

This paper is organized as follows. In Sec. [TTl we be- 
gin by presenting the methods used for determining the 
electronic structure, the geometries, and the transport 
properties of the molecular contacts. In Sec. IIIII we then 
describe the optimized geometries and the corresponding 
vibrational modes. Following this. Sec. [iVl relates in de- 
tail the transmission properties of the junctions to the 
relevant molecular orbit als. The section ends with an 
analysis of the evolution of the conductances when the 
contacts are opened. In particular, we analyze the tun- 
neling regime, which is discussed further in App. [Al In 
Sec. [Vlwe conclude with the discussion of our results. 



II. METHOD 

In this section we briefiy present our methods for com- 
puting the electronic structures, geometries, vibrational 
modes, as well as the charge transport characteristics of 
atomic contacts. A more detailed account will be given 
elsewhere 

We describe our contacts at the level of density 
functional theory (DFT) as implemented in the quan- 
tum chemistry package TuRBOMOLE v5.7 (Refs. flsHTol . 
[lOIIll). In particular, we use the "ridft" module^^^ 
with the BP86 exchange correlation functiona l^^i^^ and 
a Gaussian basis set of split valence plus polarization 
(SVP) quality;^ To be precise, polarization functions are 
present on all non-hydrogen atoms within the employed 
basis set. The total energy is converged to a precision 
of 10~^ Hartree. Geometry optimization is performed 
until the maximum norm of the Cartesian gradient has 
fallen below 10~^ atomic units. Vibrational modes in 
our contacts are determined using the "aoforce" module 

of TURBOMOLE^^i^ 

The investigated atomic contacts consist of two atom- 
ically sharp electrode tips and a single oxygen molecule 
connecting them (Fig.[Tj). The electrodes are Al fee pyra- 
mids oriented in the (111) direction with a lattice con- 
stant of 4.05 A and consist of at least four layers. While 
the outer layers are kept fixed all the time, the inner two 
layers of the freestanding pyramids are relaxed to find 
an energy minimum. Afterwards the atoms in the Al 
pyramids are kept fixed. Between two such electrodes 
we place one representative oxygen molecule (O, O2, or 



O3). The molecule is then relaxed between the fixed 
electrodes, leading to a stable configuration. Starting 
from this first configuration, we move the electrodes step- 
wise further away from each other and relax the oxygen 
molecule again in each step. For each relaxed structure 
we determine the vibrational modes of the molecule con- 
nected to the leads. 

To compute the charge transport through a molecule 
between metallic electrodes, we apply a method based on 
standard Green's function techniques and the Landauer 
formula expressed in terms of a local non-orthogonal 
basisJii^i^i^i^^i^ Due to the locality of the basis, the 
system is conveniently separated into three parts, leading 
to a single-particle Hamiltonian (or Fock matrix) of the 
molecular junction that has the following form 



H = I HcL Hcc HcR 

yHuL Hrc HrrJ 



(1) 



Here Hcc describes the central system (C) consisting of 
the molecule and the tips of the electrodes and Hxx 
with X = L^R describe the left (L) and right {R) 
electrodes. The matrices Hlc = Hcl ^^c. give the 
hopping elements between the different subsystems, but 
the hoppings between L and R are assumed to vanish 
{Hlr = H RL = 0)- This assumption has been ver- 
ified in all our calculations, where the Hamiltonian ele- 
ments connecting L and R were always smaller than 10~^ 
Hartree. The overlap matrix S of the non-orthogonal ba- 
sis has a similar structure. 

The conductances of the molecular junctions are de- 
termined by means of Green's function techniques. The 
low-temperature zero-bias conductance G = dl / dV\y=Q 
is given by the Landauer formula 



G = GoTi[t{EF)t\EF)] = GoT{Ef), 



(2) 



where t{EF) is the transmission matrix of the molecular 
junction evaluated at the Fermi energy Ef- The total 
transmission at any energy E can be written as T{E) = 
X]^i^n(^), where Tn{E) are the transmissions of the 
individual eigenchannels, defined as the eigenvalues of 
t{E)t\E). The transmission matrix can be calculated 
in terms of Green's functions of the molecular junction 
as follows^ 



t{E) =2r^/^(^)G 



cc 



{e)t'^\e). 



(3) 



Here 



G'cciE) = [EScc - Hcc - ^1{E) - (4) 

is the retarded Green's function of the central system, 
and Tx{E) = —lm[E,x{^)] scattering rate ma- 

trices. These depend on the self-energy matrices 

S^(i=;) = {Hcx-EScx)9'xxiE){Hxc-ESxc), (5) 

where g^x{E) is the retarded surface Green's function 
of the electrode X = L.R, which are modeled as ideal 
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FIG. 1: (Color online) Ali7-Oa;-Ali7: Optimized contact geometries and their tip-to-tip distances as obtained by stretching an 
aluminum (blue) atomic contact in the presence of oxygen (red) molecules. The Al pyramids consist of 7, 6, 3, and 1 atoms 
in the layers from the outside to the inside. 
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TABLE I: Distances between the aluminum tip atoms and the lengths of the individual Al-0, 0-0, and 0-Al bonds for the 
contact geometries depicted in Fig. [1] The order from top to bottom corresponds to an order from left to right in the figure. 



surfaces. These self-energies contain all the information 
about the electronic structure of the electrodes and their 
coupling to the central system. 

The electronic structure of the central system and its 
coupling to the electrodes are obtained from the calcu- 
lated geometries (Fig. [1]). Let us refer to these geome- 
tries, which consist of the two pyramids bridged by a 
molecule, as "contact geometries". These contact ge- 
ometries are divided into three parts. The central part, 
corresponding to region C in the transport calculation, 
consists of the molecule and the inner two layers of the 
pyramids. The excluded outer layers on the left and right 
sides are used for calculating the hoppings and overlaps 
between the central region and the electrodes as needed 
in Eq. 

To obtain the surface Green's functions g^^{E)^ we 
compute separately the electronic structure of a spher- 
ical Al fee cluster with 555 atoms. From this we ex- 
tract the Fock and overlap matrix elements between the 
atom in the origin of the cluster and all its neighbors 
and, using these "bulk parameters", construct a semi- 
infinite periodic crystal. The surface Green's functions 
are then calculated from this crystal using the so-called 
decimation technique. We have checked that the elec- 
trode construction in the employed non-orthogonal SVP 



basis set has converged with respect to the size of the 
Al cluster, from which we extract our parameters jii^ In 
this way we describe the whole system consistently within 
DFT, using the same basis set and exchange-correlation 
functional everywhere. 

The Fermi energy Ep of the coupled system, namely 
the C part connected to the semi-infinite electrodes, re- 
quired to determine the conductance [Eq. ([2|)], is assumed 
to be given by the Al leads. We obtain its value as the 
average of the energies of the highest occupied (—4.261 
eV) and the lowest unoccupied (—4.245 eV) orbitals of 
the bulk cluster. In this way we find Ef = —4.25 eV. 
The same quantity for the contact geometries deviates 
from this bulk value by less than 0.08 eV. Let us stress 
that the Fermi energy is only used to read off the conduc- 
tance from a transmission function T{E)^ but does not 
enter into our calculations otherwise. 



III. CONTACT GEOMETRIES AND 
VIBRATIONAL MODES 

We simulate the opening of a contact for the three most 
basic forms of oxygen (O, O2, and O3). Fig. [T] shows a 
selection of the relaxed geometries obtained, the starting 
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FIG. 2: Ali7-Occ-Ali7: Energy of the vibrational modes of the oxygen molecule connected to Al pyramids as a function of 
the distance d between the Al tip atoms. The arrows and plus/minus signs denote the directions of relative displacements 
of the atoms in a given mode: and ^ denote motion along the axis of the contact while +,— and t,| denote motion in 
the two directions perpendicular to the axis. 



configuration, one halfway, the one just before, and the 
one directly after the breaking. In the case of a single O 
atom, the O sits at first on the side of the Al dimer con- 
tact. Upon stretching this atom moves onto the contact 
axis and stays symmetrically bonded to both pyramids 
until the contact breaks [Fig. |l(a)] . When the contact is 
broken, the oxygen atom binds to one electrode with a 
distance of approximately 1.7 A on top of the pyramid 
and stays bonded like this. However, this binding posi- 
tion is only metastable. If one distorts the position of 
the oxygen atom slightly, e.g. places it at some distance 
from the contact axis and optimizes the contact geome- 
try again, the atom goes to one side of the pyramid in 
the same way as for O2 and O3 described below. This 
position with the oxygen bonded to one side of the pyra- 
mid is energetically favored by approximately 2.5 eV. In 
the case of O2, the molecular axis of the oxygen molecule 
is in the beginning almost perpendicular to the contact 
axis and rotates afterwards into this axis [Fig. [1(b)] . For 
O3, the molecule is twisted around the contact axis. This 
twist contracts towards the contact axis when the contact 
is stretched [Fig. l(c)| . In the cases of AI-O2/O3-AI, the 
O atoms move always to the sides of the pyramids when 
the contact is broken. The evolution of the distances be- 
tween the atoms while opening the contact is shown in 
Tab.m 

As mentioned above, in order to compare with results 
from inelastic transport experiments, we also calculate 
the vibrational modes of the oxygen molecules in the con- 
tact. Since the atomic mass of Al is almost twice as large 
as that of O, it is a reasonable first approximation to 
neglect the motion of the Al atoms, which is what our 
model assumes. The energies of the vibrational modes 
for the different structures are shown in Fig. [2l For these 
structures, the modes can be classified roughly as being 
longitudinal or transverse, according to the motions of 
the atoms relative to the axis of the contact. The mode 
energies are in the range of 50 meV to 140 meV for Al-0- 
Al. The two transverse modes are more or less constant 



around 50 meV, while the longitudinal mode increases 
from 120 meV to a maximum of 140 meV at a distance 
of 3.3 A and drops afterwards almost linearly. For AI-O2- 
Al, the mode energies are in the range of 40 meV to 150 
meV. The two lowest modes are transverse and constant 
in energy. The highest mode is longitudinal and increases 
from 135 meV to 150 meV at 4.5 A and drops afterwards 
to 120 meV. The other three modes lie in the range of 
70 meV to 110 meV. For AI-O3-AI, the lowest (trans- 
verse) modes lie constant at an energy of approximately 
25 meV. The highest (longitudinal) mode increases first 
to its maximum of 160 meV at 5.4 A and decreases af- 
terwards. Another mode lies close to this highest one. 
All other modes have energies between 50 meV and 120 
meV. 

The behavior of the vibrational energies can be under- 
stood by re-examining the atomic distances in the geome- 
tries (Tab. [J). Despite the fact that the atomic contacts 
are being stretched, it happens for example that in the 
beginning some atoms move closer together. Therefore 
certain bonds strengthen due to this bond shortening, 
and corresponding modes may show an initial increase in 
their vibrational frequencies. 



IV. TRANSPORT PROPERTIES AND 
ELECTRONIC STRUCTURE 

Before we determine the transport properties for the 
different forms of oxygen between aluminum electrodes 
we will show the results for a pure Al atomic contact. 
The transmission as a function of energy and the corre- 
sponding contact geometry are plotted in Fig. [3l The 
conductance of this contact is close to IGq and is carried 
by 3 channels, as expected i^^^^ If oxygen sits be- 
tween the two electrodes, the transmission changes signif- 
icantly. The resulting T{E) curves for the contact geome- 
tries of Fig. [1] are shown in Fig. [H We note that when- 
ever we speak about numbers of channels in this work, we 
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FIG. 3: (Color online) Al-Al: Transmission as a function of 
energy. Both the total transmission (T) as well as its chan- 
nel contributions (Tn) are shown. The vertical dashed line 
indicates the Fermi energy of bulk Al. The corresponding 
contact geometry with a distance between the Al tip atoms 
of 3.4 A is depicted to the right. 



mean channels that contribute at least 5% to the total 
conductance. This criterion approximately agrees with 
what can be resolved experimentally using superconduct- 
ing electrodes. 

The total transmissions T{E) = ^^Tn{E) and the 
transmission eigenvalues Tn{E) of the individual channels 
of Al-O-Al [Fig. [T(a)] are plotted as a function of energy 
E in the left column of Fig. [H The peaks labeled 1- 
3 below the Fermi energy are due to oxygen. This can 
be seen from the local density of states (LDOS) shown 
for a tip-to-tip distance of 3.7 A in Fig. [5] (Ref. [iO) as 
well as from the orbitals at the peaks shown in Fig. [6l 
Peak number 1 corresponds to the Pz orbital of O {z is 
the transport direction). Peaks 2 and 3 become almost 
degenerate as soon as the O is on the contact axis. These 
peaks are due to the Px and Py orbitals of O. Peak 1 
is broader and shifted to lower energies than peaks 2 
and 3 due to the stronger binding of the Pz orbital to 
the Al states. The peaks above the Fermi energy are 
due to aluminum. The number of channels contributing 
to the conductance is at least three. When the contact 
is opened, peaks 1-3 move, up to a tip-to-tip distance 
of 3.3 A, away from and afterwards towards the Fermi 
energy. This behavior is caused by the length of the Al- 
O bond. 



For AI-O2-AI [Fig. 1(b) , the total transmission and the 
channel transmissions are shown in the central column of 
Fig. m The peaks in T{E) are caused by the oxygen 
as can be seen from the LDOS in Fig. [5] for a tip-to-tip 
distance of 4.5 A. Peaks number 4 and 5 originate from 
the anti-bonding tt orbitals and peak 6 from the anti- 
bonding a orbital (Fig. [6]). The number of channels at 
the Fermi energy is two for the intermediate tip-to-tip 
distances between 4.6 A and 4.9 A, but otherwise three. 
Peaks 4 and 5 shift upwards in energy without crossing 



the Fermi energy while the contact is stretched. Peak 6 
moves first away from and then, after a tip-to-tip distance 
of 4.5 A, towards the Fermi energy. The movement of this 
peak can be understood in terms of the behavior of the 
distance between the O atoms, while the length of the 
Al-0 bond is mainly responsible for the displacement of 
peaks 4 and 5. 

For AI-O3-AI [Fig. l(c)| , the transmissions are plotted 
in the right column of Fig. [H The peaks shown have 
their origin in the ozone (Fig. [5] and [6]). Peaks 8 and 9 
move away from the Fermi energy with stretching up to 
a tip-to-tip distance of 5.4 A and approach afterwards 
the Fermi energy. Peak 7 shifts continuously towards the 
Fermi energy. The movement of peaks 8 and 9 again 
reflects the distances between the oxygen atoms, and the 
shift of peak 7 changes in the Al-0 bond length. The 
current is carried, up to a tip-to-tip distance of 5.4 A, by 
three channels and afterwards by two channels. 

In Fig. [Tithe conductance is plotted as a function of the 
tip-to-tip distance for pure Al and for the three different 
oxygen molecules in Al contacts. For pure Al contacts we 
obtain a conductance close to IGq with a positive slope 
shortly before rupture of the contact, in agreement with 
previous investigations of Al atomic contacts J^i^i^i^ 
In the presence of oxygen, the conductances show the 
same qualitative features for the three molecules (O, O2, 
and O3). The conductance first decreases with stretch- 
ing up to a certain point. After this point it increases 
rapidly until the contact breaks and the conductance 
drops. This trend can be understood by re-examining the 
evolution of the energy dependence of the transmission 
upon stretching (Fig. [4|). Thus for instance, the increase 
of the conductance originates from the fact that the high- 
est occupied molecular orbitals move close to the Fermi 
energy. Quantitatively the conductance depends on the 
number of oxygen atoms in the molecule. The more oxy- 
gen atoms, the lower the conductance. In all cases, the 
transport is influenced by a charge transfer from the Al 
to the oxygen molecule (around — 0.6e for all the three 
molecules, as determined by a Mulliken population anal- 
ysis) leading to an occupation of molecular orbitals which 
would be empty in an isolated oxygen molecule. There- 
fore, current is carried by formerly unoccupied p orbitals. 

We also investigate the conductance in the tunneling 
regime. With this term we refer to the region of large 
tip-to-tip distances where the conductance exhibits an 
exponential decay G{d)/Go cx e~^^ with an inverse decay 
length /3. In experiments there might be any number of 
O molecules covering the surfaces of the Al pyramids, but 
we calculate representatively the tunneling conductance 
of AI-O2-AI (not shown in Fig. [7j), where the O atoms 
sit separated on both pyramids [see the rightmost part of 
Fig. [1(b)]. We compare this to the tunneling conductance 
of pure Al contacts (shown in Fig. [7j). For pure Al, the 
inverse decay length turns out to be /3 = 2.1 A~^ and for 
Al with an O atom on both pyramids (3 = 2.3 A~^. Thus 
the decay of the conductance is somewhat faster when 
oxygen is present. 
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FIG. 4: (Color online) Al-Ox-Al: Transmission as a function of energy for the structures shown in Fig. [T] The tip-to-tip 
distance increases from the top to the bottom and is given in the panels. The vertical dashed lines indicate the Fermi 
energy. 



Physically, the value of (3 is determined by the shape 
of the potential barrier between the tips, especially its 
apparent height with respect to Ep (Ref. 43). The bar- 
rier depends sensitively on the electronic structure of the 
junction, including the charge-transfer effects between 
the different atoms close to the tips. In terms of the 
local-basis picture, all this information is contained in 
the spectral functions of the tips and the hopping matrix 
elements connecting them. The mathematical details of 
this interpretation are discussed further in App. [Xj In 
this picture, an important role is played by the radial 



decay properties of the orbitals of the tip atoms that are 
relevant for transport. Roughly speaking, for O these 
orbitals are the 2p ones, while for Al they are of the 35 
and 3p type. The former orbitals decay faster and thus 
it is reasonable that inverse decay length j3 is larger in 
the presence of oxygen. Although the oxygen-induced 
change in /3 in this particular example is not very large, 
we suggest that, in general, changes in the decay rate 
of the tunneling conductance may provide another possi- 
bility to check experimentally whether gas molecules are 
present in the atomic contact or not. 
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FIG. 5: (Color online) Al-Ox-Al: LDOS corresponding to the second row in Fig. |4l The solid (black) curve is for the 
region including the Al tip atoms and the O atoms, the dash-dotted (red) one is for one of the oxygen atoms, and the 
dashed (blue) curve is for one of the Al tip atoms. The vertical dashed (orange) line indicates the Fermi energy. 




FIG. 6: (Color online) Alir-Ox-Alir: Orbitals of the contact geometries of Fig. [T] The orbitals 1-3 belong to 
AI17-O-AI17 with a tip-to-tip distance oi d — 3.7 A, the orbitals 4-6 to AI17-O2-AI17 with d — 4.5 A, and the 
orbitals 7-9 to AI17-O3-AI17 with d = 5.6 A. The energies of the orbitals are -10.537 eV (1), -7.238 eV (2), -7.176 
eV (3), -6.274 eV (4), -5.606 eV (5), -2.133 eV (6), -5.197 eV (7), -2.175 eV (8), and +0.153 eV (9), where the 
number of the orbital is given in brackets. 



V. DISCUSSION AND CONCLUSIONS 

Let us point out that we have checked that the results 
presented in this work are robust with respect to the size 
of the contacts. In particular, the charge transport char- 
acteristics do not change significantly when the size of 
the pyramids in the contact geometries is increased, nor 
do they depend sensitively on the partitioning into the L, 
C, and R regions. We illustrate this fact in Fig. [3 where 
we compare the transmission functions T{E) of the Al- 
O2-AI contact of Fig. 1(b) | at a tip-to-tip distance of 4.5 
A. We vary the numbers of layers in the Al pyramids 
and divide the system differently into the three regions. 
Clearly, the essential features of the results remain unaf- 
fected by these modifications. Such a robustness is ob- 
served as long as the assumption about vanishing matrix 
elements between the L and R regions is fulfilled, and 
there are enough layers (at least two in the case of Al) 
left in the L and R parts of the pyramids for coupling 
the C part to the surfaces. 



In order to investigate structures where the effect of 
O molecules on the transport properties is at its largest, 
we have placed the oxygen molecules initially between 
the pyramids and kept the Al electrodes fixed in the 
geometry-optimization process. On the other hand, we 
know that relaxing parts of the pyramids together with 
the molecule can lead to geometrical structures differing 
considerably from the ideal contacts shown in Fig. [TJ An 
example of such a geometry is depicted in Fig. [9l More- 
over, the formed structures depend on the initial place- 
ments of the atoms close to the tips and on the num- 
ber of oxygen molecules present. Different realizations 
of atomic configurations can also be observed in exper- 
iments, where various kinds of I-V characteristics have 
been measured. However, as the transmission function 
in Fig. [9] shows, even more complex geometries can ex- 
hibit results whose main features resemble those of the 
simplified structures discussed above. The transmission 
function shown in this figure should be compared with 
the ones for the Al-O-Al contact in Fig. [H In particu- 
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Distance d (A) 



FIG. 7: Al-Occ-Al: Conductance as a function of the 
distance between the Al tip atoms. 




Energy E (eV) 



FIG. 8: (Color onhne) Robustness of the transmission 
for an AI-O2-AI contact at a tip-to-tip distance of 4.5 A 
[Fig. mb)]. Shown are T{E) curves for an increased 
number of Al atoms in the pyramids as well as for dif- 
ferent partitionings of the contact into L, C, and R 
regions. The solid black line is the transmission of 7- 
6-|3-l-02-l-3|-6-7 (Fig.[ll). The numbers stand for the 
Al atoms in each layer of the pyramid and "|" indicates 
how the system is divided into the three different re- 
gions. The dash-dot-dotted Une corresponds to the sys- 
tem 6-12-|10-6-3-l-O2-l-3-6-10|-12-6, the dashed line to 
6-12-10-|6-3-l-O2-l-3-6|-10-12-6, and the dash-dotted 
line to 6-12-10-6-|3-l-O2-l-3|-6-10-12-6. The vertical 
dashed line indicates the Fermi energy. 



lar, similarly to those results, there are three peaks (1-3) 
below the Fermi energy and the conductance is around 
O.2G0. 

A similar approach to investigate oxygen in Al junc- 
tions was chosen by Jelinek et al. (Ref. 1421 ). Those au- 
thors found structures where an oxygen atom is incor- 
porated into the electrodes and their electrodes are con- 




Energy E (eV) 



FIG. 9: (Color online) AI-O2-AI: Example of a pos- 
sible geometrical configuration when the Al tip atoms 
are included in the relaxation process in addition to 
the molecule, and the corresponding transmission as a 
function of energy. The vertical dashed line indicates 
the Fermi energy. 



nected by an aluminum dimer in the final stage before 
the contact breaks. The conductance in their struc- 
tures, therefore, is dominated by the aluminum dimer 
and the influence of the oxygen atom on the conduc- 
tance is weak. We nevertheless believe that structures 
with oxygen molecules bridging the gap between Al elec- 
trodes can be found in experiments. A situation corre- 
sponding to our simulations would most likely arise when 
the oxygen is introduced into open contacts that are sub- 
sequently closed and then reopened. The experimental 
results can then be compared to our predictions. 

In conclusion, we studied the charge transport prop- 
erties of aluminum atomic contacts in the presence of 
oxygen molecules (O, O2, and O3). We obtained stable 
geometrical structures for the molecules in the junctions, 
and determined the molecules' vibrational modes in the 
presence of the aluminum electrodes. In principle these 
modes should be measurable with the methods of point 
contact spectroscopy. Moreover, we analyzed the evo- 
lution of the conductance for the AI-O/O2/O3-AI con- 
tacts in stretching processes, mimicking a conductance 
curve measured in experiments while opening a break 
junction. These results show a typical behavior for all 
the different molecules studied. First the conductance 
decreases, until it starts to increase shortly before the 
contact breaks. The value of the conductance was found 
to decrease with the number of oxygen atoms. Our ob- 
servations can be understood by the interplay between 
mechanical and electronic properties. Changes in the 
distances between the atoms result in modifications of 
the electronic structure. The shifts in energy positions 
of individual molecular orbitals and their corresponding 
peaks in the transmission function can in turn explain the 
behavior of the conductance in the stretching processes. 
Finally, we studied the decay lengths of the conductance 
in the tunneling regime, and found the decay to be faster 
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in the presence of oxygen than for pure aluminum con- 
tacts. This effect may provide a way for detecting gas 
molecules in metallic atomic contacts. 
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APPENDIX A: CONDUCTANCE FORMULA IN 
THE TUNNELING LIMIT 

In the actual conductance calculations presented in 
this paper, we always use the formula given in Eq. ([2]). 
For a better understanding of the tunneling regime, it 
is nevertheless useful to cast Eq. (j2j) into another well- 
known form. Let us assume that the C part of the system 
can be divided into two regions 1 and 2, where region 
1 (2) is not coupled to the R (L) lead through direct 
hoppings or overlaps. Furthermore, regions 1 and 2 are 
connected to each other by ti2 = H12 —ES12. Equation 
([2]) can then be written in the form 

G = GoTr[G^2(rfl)22G'^i(rL)ii], (Al) 

where G21 = (^12)^5 energy-dependent quan- 

tities are evaluated at Ep. Let us denote by g^i and 
the Green's functions for the regions 1 and 2 in the 
absence of ti2. For example gii{E) = [ESu — Hu — 



(5]2)ii] ^. We also define the unperturbed spectral den- 
sity matrix An = i{g{^ -9ii), where g^^ = (c/n)^ with 
a similar definition for A22. We note that Hn and H22 
are simply cut out from the full, self-consistent Hoc and 
thus still contain indirect effects of ti2. 

The full propagator between 1 and 2 is given by the 
Dyson equation"^^ = g\iti2Gl2 = 9uTi2g22- 

the last stage, we defined the t matrix T12 due to the 
perturbation ti2 (Ref. [46l ). Inserting these into Eq. 
()Aip and using the identities gii{'i^L)iigii = ^11 and 
922{^R)'22g22 — ^22, wc find the following formula^'^ 

G = GoTr[Ti2A22T2iAn]. (A2) 

Far in the tunneling regime, the quantities t\2 are small, 
such that the lowest-order approximation T\2 = ti2 
should be valid. Furthermore, to lowest order, the elec- 
tronic structures of the tips (as represented by Hu and 
^22, and hence An and A22) also become independent 
of the tip-to-tip distance. 

The decay rate of the conductance with increasing dis- 
tance between the tips {P) is now seen to be determined 
by two factors. The first is due to the decay of the cou- 
plings ti2. If one uses the Wolfsberg-Helmholz (or "ex- 
tended Hiickel") approximation, then H12 and hence ti2 
are related in a simple way to the overlap S12 (Ref. 1481 ). 
This emphasizes that the major effect to the decay of 
ti2 is simply due to the associated local basis functions, 
which decay with varying rates. The second factor is 
due to the spectral functions An and A22. They weight 
the various components of ti2 differently, and thus fi- 
nally determine the decay rate /3. The physical quantity 
controlling this rate is still the effective tunneling bar- 
rier, in particular its height with respect to Ep, because 
the latter determines the barrier-penetration lengths of 
the electronic eigenstates of the tips. This information 
is contained in the matrices An and A22 in the form of 
the weights of the different local basis functions of the 
tip atoms. 
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